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ABSTRACT 

We present new multi-band transit photometry of three small {Rp ^ 6 f?®), short- 
period (P 6 days) Kepler planet candidates acquired with th e Gran Telescopio Ca- 
nari as. These ob s ervat ions supplement the results presented in IColon fc Ford! (j201ll) 
and IColon et al.l l|2012f) . where we used multicolor transit photometry of five Kepler 
planet candidates to search for wavelength-dependent transit depths and either vali¬ 
date planet candidates or identify eclipsing binary false positives within our sample. 
In those previous studies, we provided evidence that three targets were false positives 
and two targets were planets. Here, we present observations that provide evidence 
supporting a planetary nature for KOI 439.01 and KOI 732.01, and we find that KOI 
531.01, a 6 P® planet candidate around an M dwarf, is likely a false positive. We also 
present a discussion of the purported “sub-Jovian desert” in the orbital period-planet 
radius plane, which cannot be easily explained by observational bias. Both KOI 439.01 
and KOI 732.01 are likely planets located within the so-called desert and should be 
investigated with further follow-up observations. As only ^30 of the '^3600 currently 
active Kepler planet candidates are located within the sub-Jovian desert, it will be 
interesting to see if these candidates also survive the vetting process and fill in the 
gap in the period-radius plane. Confirming planets in this regime will be important for 
understanding planetary migration and evolution processes, and we urge additional 
follow-up observations of these planet candidates to confirm their nature. 
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1 INTRODUCTION 


With the discovery of nearly 1000 planets and several thou¬ 
sand additional planet candidates awaiting confirmation Q 
the Kepler space mission has set the stage for in-depth 
studies of different planet populations. In order to study 
these populations, the r eliability of the Kep l er sample must 
be tak e n into accoun t . iMorton fc JohnsonI (l20 1lh. iMortorJ 
ll2012ll . iFressin et al.l (l2013l l. and Desert et al.l Eolsi l find 
a low global false positive rate for Kepler candidates of 
~10%, suggesting tha t most of the remainin g candidates are 
real planets. However. ICouehlin et al.l (l2014h estimate a 35% 
false positive rate f or Kepler candi d ates d ue to various types 
of con tamination. ISanterne et al.l ll2012h and IColon et al.l 
ll2012ll also suggest higher false positive rates of 35% and 
~50—60% for specific subsets of Kepler candidates, such as 
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small candidates in single syste ms with orbital pe riods ap¬ 
proximately less than six days llColon et al.l 12012 1. On the 
other hand, nearly all can didates in multi-candidate systems 
are e xpected to be real llLissauer et al. I I2OI4I : I Rowe etld] 
l2014h . 

Given this range of false positive rates, it is unclear 
how significantly false positives affect statist ical studies of 
differ ent subsets of Kepler candidates (e.g. [Howard et al.l 
l2012li . In particular, there is an apparent excess of single¬ 
candidate systems (relative to multi-candidate systems), but 
it remains to be seen what fraction of the single-candidate 
systems are real. For this reason, candidates believed to be in 
single systems are the focus of this paper. Such candidates 
with short-periods are particularly interesting, given that 
some are located within an apparent sub-Jovian desert in 
the orbital period-planet radius plane (ISzabo fc KissI 1201 ll : 
iBeauee fc Nesvorn^l2013l i. 

In IColon fc FordI (1201 ll i and IColon et al.l ll2012li (here¬ 
after, Paper I and Paper H), we presented multi-wavelength 
transit observations of five small, short-period Kepler planet 
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Table 1. KOI Star Properties 


KOI 

KIC 

Kp 

Te// (K) 

b (deg) 

439 

12470954 

14.313 

5438 

13.06 

531 

10395543 

14.418 

4030 

16.61 

732 

10265898 

15.342 

5379 

16.15 


All values are from the NASA Exoplanet Archive 
(http://exoplanetarchive.ipac.caltech.edu/ accessed on 2015 
April 4). The KIC number is the Kepler Input Catalog number 
for each target. Here, b is the Galactic latitude of the KOI host 
star. 

candidates (KOI 565.01 and KOI 225.01, 420.01, 526.01, and 
1187.01) acquired with the 10.4-m Gran Telescopic Canarias 
(GTC). A complete description of our target selection crite¬ 
ria is given in Paper II and is summarized here, again with 
only single-candidate systems being considered: 

• orbital period (P) ^ 6 days 

• planet radius (Rp) 6 7?® 

• transit depth (5) ^ 500 ppm 

• transit duration (r) ^ 2.5 h 

• Kepler magnitude (Kp) ^ 15.5 

Gurrently, there are 56 active KOIs that satisfy these crite¬ 
ria, including KOI 420.01 and 526.01 (the two systems we 
validated as planets in Paper II). In these previous papers, 
we looked for wavelength-dependent transit depths for each 
target, which is suggestive of an eclipsing binary false posi¬ 
tive scenario assuming the target is not a system composed 
of two (or more) equal mass or equal temperature eclipsing 
stars (which would reveal no color difference). Of hve tar¬ 
gets observed, we identihed three as having achromatic light 
curves and validated two as real planets. The resulting high 
false positive rate for our sample was not surprising, given 
the large percentage of eclipsing binaries with short orbital 
periods that can mimic planetary transits (particularly when 
diluted by other nearby or unresolved stars). We note that 
we did not inspect the Kepler light curves of our targets 
in detail prior to observing them. We assumed the targets 
had been appropriately vetted prior to being added to the 
candidate list, so we did not search for secondary eclipses 
in the Kepler data nor did we identify candidates with par¬ 
ticularly V-shaped light curves. However, targets displaying 
such attributes would be more likely to be false positives 
than a typical candidate, so we encourage similar studies to 
perform detailed inspections of Kepler light curves prior to 
collecting any follow-up data. 

In this paper, we present multi-wavelength transit ob¬ 
servations from the GTG of three additional small, short- 
period Kepler planet candidates (KOI 439.01, 531.01, and 
732.01). Tables[T]and[2]list some relevant properties for each 
of these targets. In Sections [5] and |3] we present a summary 
of our observations and our data reduction and analysis. 
Results for each target are presented in Section [I] A dis¬ 
cussion of our results for each target, as well as an updated 
investigation of trends in the false positive rate with differ¬ 
ent planetary and stellar parameters, is given in Section [5] 
We also include a discussion of candidates located within 
the apparent sub-Jovian desert in Section [5] and we offer 
concluding remarks in Section |6] 


2 OBSERVATIONS 

We used OSIRIS on the GTC to acquire near-simultaneous 
multicolor transit photometry of three Kepler planet candi¬ 
dates: KOI 439.01, KOI 531.01, and KOI 732.01. This facility 
is described in further detail in Papers I and II. As in Paper 
n, we alternated between two broadband order sorter filters 
during each transit observation: 666 ± 36 nm and 858 ± 
58 nm. All observations were conducted in queue (service) 
mode and used 1x1 binning and a fast readout rate of 500 
kHz. We describe specific details for each observation in the 
following sections. 

2.1 KOI 439.01 

The 2012 June 11 UT transit of KOI 439.01 was observed un¬ 
der photometric conditions during gray time. Observations 
began at 00:18 UT and ended at 04:30 UT on 2012 June 11, 
during which time the airmass varied from 1.08 to 1.43. At 
01:40 UT there was a technical problem with the primary 
mirror, and the time series was interrupted. Observations 
resumed at 02:04 UT. The seeing was extremely variable in 
the first part of the night, ranging from 0.9 to 1.7 arcsec 
(7.1 to 13.4 pixels). A slight defocus was implemented after 
03:00 UT due to an overall improvement in the seeing. The 
integration time was set to 60 s for both filters, and a sin¬ 
gle window of 1900 x 2100 pixels located on one CCD chip 
was read out. The resulting dead time between each expo¬ 
sure was 21 s. Autoguiding kept the target centroid stable 
to within 7 pixels. 

2.2 KOI 531.01 

We observed the 2012 September 5 ut transit of KOI 531.01 
in photometric/clear conditions during bright time. The air- 
mass ranged from 1.07 to 1.53 during the observations, 
which took place from 22:01 UT (2012 September 5) to 01:20 
UT (2012 September 6). The seeing was unstable (varying 
between 0.8 and 1.5 arcsec or 6.3 and 11.8 pixels), so the fo¬ 
cus was adjusted several times to avoid saturating the target 
star. Still, the target was saturated in four images, and we 
exclude these from our analysis. The integration time used 
for both filters was 23 s. The full CCD chip where the target 
was located was read out, resulting in a dead time between 
exposures of 29.5 s. The target centroid drifted by < 4 pixels 
during the observations. 

2.3 KOI 732.01 

Observations of the 2012 May 27 UT transit of KOI 732.01 
took place under clear conditions and during gray time. The 
observations began at 00:42 UT and ended at 05:12 UT on 
2012 May 27. In that time the airmass ranged from 1.05 to 
1.37 and the seeing varied between 0.8 and 1.0 arcsec (6.3 
and 7.9 pixels). No defocus was applied. An integration time 
of 60 s was used for both filters. Both CCD chips were read 
out, but in our analysis we only use comparison stars that 
are located on the same chip as the target to minimize any 
systematic differences in the properties of the two chips. The 
corresponding dead time between exposures was 34 s. Auto¬ 
guiding kept the target centroid stable to within 5 pixels, but 


© 2015 RAS, MNRAS OOO.ITlfTOl 





Vetting Kepler planet candidates 3 


Table 2. KOI Candidate Properties 


KOI 

P (days) 

T (d) 

b 

R-p/Ri, 


<5 (ppm) 

EB Catalog 

439.01 

531.01 

732.01 

1.90221±4.1e-07 

3.68747±3.9e-07 

1.26026±6.6e-07 

0.09199±2.875e-04 

0.04198±7.500e-04 

0.07809±6.542e-04 

0 0433+l-69e — 01 
u.U‘iOO_4 32e_02 

n — 92 

U.yu/ _i 09e_02 
n — 92 

0.0435l3-04e-04 

0.0868l®;«:i“3 

0.035313-80-05 

Q Qo-|-1.19e—00 
'^•^^-2.50e-01 
5.69±6.60e-01 

i:?Q+l-30e—00 

'J-'J^_2.90e-01 

2341±7.900e-|-00 

5295±2.550e+01 

1207±1.030e+01 

yes 

yes 

no 


All values are from the NASA Exoplanet Archive (http://exoplanetarchive.ipac.caltech.edu/ , accessed on 2015 April 4). Here, r is the 
transit duration, b is the impact parameter, and 5 is the transit depth as measu red in the Kepler ban dpass. The last column indicates 
if the KOI is listed in the third version of the Kepler eclipsing binary catalog llSlawson et al.ir201ll: http://keplerebs.villanova.edu/; 
accessed on 2015 April 4). 


there are a few outliers towards the end of the observations 
where the centroid drifted by as much as 11 pixels. 


3 DATA REDUCTION AND ANALYSIS 

All images were reduced using software written in GDL0 
We performed bias subtraction and flat fielding on each sci¬ 
ence image prior to performing aperture photometry. Dome 
flats acquired for each filter were used, and we corrected 
for non-uniform illumination by the dome lamp by remov¬ 
ing the large-scale illumination pattern in the final co-added 
flat field for each bandpass through smoothing. We note that 
no flats were taken for the 858/58 filter for the KOI 531 ob¬ 
servations. We therefore use the 858/58 fiats from the KOI 
732 observations in our analysis of the KOI 531 data, as 
both observations utilized the full CCD chip. Although these 
two observations were taken about three and a half months 
apart, we do not believe our analysis was significantly af¬ 
fected by using “old” flats on the KOI 531 858/58 data. In 
particular, we find that the standard deviation of the time- 
binned residuals for the KOI 531 light curve at 858-nm is 
consistent with the trend expected for white Gaussian noise 
(despite systematics that are clearly seen in the light curves; 
0 . 

Circular aperture photometry was performed on each 
target, stars within ~20 arcsec of the target (i.e. nearby 
stars that could be the source of the transit signal), and a 
sample of comparison stars. We initially used apertures with 
a range of radii, but ultimately the aperture that yielded the 
best photometry (i.e. the lowest rms scatter outside of the 
transit) for each target was used. For KOI 439, an aperture 
of 24 and 26 pixels (3.0 and 3.3 arcsec) was used for the 
666 -nm and 858-nm data, respectively, and an annulus with 
an inner radius of 60 pixels and outer radius of 65 pixels 
(7.6 and 8.2 arcsec) was used for sky subtraction. For KOI 
531, an aperture of 26 and 24 pixels (3.3 and 3.0 arcsec) 
was used for the 666-nm and 858-nm data, respectively, and 
an annulus with inner and outer radii of 50 and 55 pixels 
(6.3 and 7.0 arcsec) was used for sky subtraction. For KOI 
732, an aperture of 25 and 34 pixels (3.2 and 4.3 arcsec) 
was used for the 666-nm and 858-nm data, respectively, and 
an annulus with inner and outer radii of 70 and 75 pixels 
(8.9 and 9.5 arcsec) was used for sky subtraction. To avoid 
including the target flux in the sky annulus of the nearby 


^ GNU Data Language; http://gnudatalanguage.sourceforge.net/ 


stars we checked as potential sources of the transit signal, 
we sometimes applied a larger sky annulus for those stars. 

We generated a light curve for each target by dividing 
the total flux measured from the target star by the total 
weighted flux of an ensemble of comparison stars. Compar¬ 
ison stars that had a similar brightness to the target and 
that did not display variability were included in the ensem¬ 
ble. For consistency, for a given target the same ensemble 
was used to generate the 666- and 858-nm light curves. For 
KOI 439, 531, and 732, ensembles of three comparison stars 
were ultimately used. We excluded any data where the peak 
counts in the target were > 50000 to avoid the saturation 
limit and stay in the linear regime of the CCD. The light 
curves were then normalized by dividing by the median flux 
ratio measured in the out-of-transit data. 

The observation epochs at mid-exposure were extracted 
from the FITS header for each image. The extracted 
UTC times were converted to Barycentric Julian Dates in 
Barycentric Dynamical Time (BJD_TDB) via an online ap- 
pleCl for consistency with the time coordinate system that 
the Kepler mission uses. 

The out-of-transit light curve for each target was then 
regressed against airmass, the target centroid position, and 
the peak counts in the target (per pixel), and a linear trend 
was removed from the full light curve for each target. For 
KOI 531, we also linearly detrended the light curve against 
the sky background and removed a quadratic trend that was 
apparent in the baseline (out-of-transit) data. 

The median photometric error for the target light curve, 
which was computed from the photon noise of the target and 
comparison stars, the noise from the sky background, and 
scintillation, is 577 and 522 ppm for the 666- and 858-nm 
light curves for KOI 439, 838 and 682 ppm for the 666- and 
858-nm light curves for KOI 531, and 693 and 613 ppm for 
the 666- and 858-nm light curves for KOI 732. In all cases, 
the photon noise of the target and reference star ensem¬ 
bles dominates the error budget. In most cases, the stan¬ 
dard deviation of the light curve residuals is larger than 
the computed photometric error, suggesting that our com¬ 
puted errors may be underestimated. We therefore perform 
a prayer-bead analysis to better characterize the uncertain¬ 
ties on each htted parameter in our light curve model. We 
discuss this further below. 

We followed a similar procedure as in Paper II to ht a 
synthetic limb-darkened light curve model lIMandel fc Agoll 


^ http://astroutils.as tronomy.ohio-state.edu/time/utc2bjd.html; 
lEastrnan et al. 111^3) 
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l2002h to the detrended data for each target. The process is 
summarized here; For each target, we fitted models to each 
light curve separately and then corrected the data against 
the best-fitting (linear) baseline slope. Light curve residu¬ 
als were computed by subtracting the model from the data 
and were used to identify and discard outlying data points 
(i.e. points greater than 3 ct from the median value of the 
residual light curve). We then fitted models to the corrected 
light curves simultaneously and forced the mid-transit time 
(tc), transit duration (from first to fourth contact; r), im¬ 
pact parameter (&), baseline flux ratio, and baseline slope 
to the same value for both light curves. The planet-star 
radius ratio (p = Rp/Ri,) and limb-darkening coefficients 
were determined separately for each wavelength. The fitting 
process involved iterating over a range of initial guesses for 
each parameter. In each iteration, each limb-darkening coef¬ 
ficient was held fixed at a specific value for each wavelength 
rather than being allowed to vary freely like the other model 
parameters. However, each iteration used a different set of 
values for the limb-darkening coefficients to account for un¬ 
certainties in the coefficients. A prayer-bead analysis was 
also performed here, where we fitted models to synthetic 
data generated from circularly shifting the residuals of each 
light curve. In Section U] we consider the distribution of the 
best-fitting parameters to the real and synthetic data when 
discussing the uncertainties for each fitted parameter in or¬ 
der to account for any sources of systematic noise. 

Finally, we computed the false positive probability 
(FPP) for each KOI presented in this paper using vespa. 
This tool is based on the Morton (2012) procedure for cal¬ 
culating FPPs for transiting planet candidates discovered by 
the Kepler mission and was recently made available for gen¬ 
eral use by the community^ The vespa tool computes FPPs 
by simulating three different astrophysical scenarios that 
could mimic a planetary transit: (1) an undiluted eclipsing 
binary (EB), a hierarchical eclipsing binary system (HEB), 
or a background (or foreground) eclipsing binary blended 
with the target (BEB). We specifically used vespa to calcu¬ 
late FPPs for each of our KOIs from the Kepler photometry 
(downloaded from the NASA Exoplanet Archive), inputting 
the right ascension and declination, orbital period, planet- 
star radius ratio (measured from Kepler', Table [2|, and the 
Kp, V, griz, and JHK magnitudes for each KOI. Given the 
lower photometric precisions achieved in the GTC observa¬ 
tions presented here, we only compute FPPs based on the 
Kepler photometry. The results of our analysis are presented 
in the following sections. 


4 RESULTS 

In Figs[T]-[3]we present the light curves and the best-fitting 
models for each target, along with the light curve residu¬ 
als and the color (666—858 nm). To compute the color we 
binned each light curve in 10-minute intervals and then com¬ 
puted the color as -2.5 log{Fem/Fs 5 s)- Here, Fme is the flux 
ratio measured in the 666 nm bandpass and Fssa is that mea¬ 
sured in the 858 nm bandpass. Thus in Figs[T]-[31 a positive 


^ http://github.com/timothydmorton/vespa 


color indicates a ‘red’ transit and a negative color indicates 
a ‘blue’ transit. 

While not presented here, we also examined the light 
curves for stars within 20 arcsec of each target to confirm 
that the target is the source of the transit signal. We found 
no transit-like signals in any nearby stars, confirming that 
in each case the target is most likely the variable star. How¬ 
ever, we cannot rule out that the signal is due to an object 
transiting an unresolved star (specifically within 2.5 arcsec 
of the target for KOI 439.01 and KOI 531.01; 1)5. H and 1)5.21) . 

In Table[3]we present the best-fitting model parameters 
for each target. The best-fitting parameters are identified as 
the median best-fitting value from the prayer-bead chain. 
The la uncertainties for each parameter are also presented 
in Table O which are computed as the standard deviation of 
each parameter fitted in the prayer-bead analysis. We also 
present in Table [4] the median ratio of the planet-star radius 
ratios measured at 666-nm and 858-nm (i.e. pme/psss)- 

In Table Owe present the results from the FPP calcu¬ 
lations. The table contains the relative probability that the 
measured light curve is either a planetary transit, an EB, a 
HEB, or a BEB. The final column contains the FPP. 

Results for each individual target are summarized here 
and are discussed in further detail in the following section. 

KOI 439.01'. We find that the transit depths in the two 
GTG bandpasses are consistent 0 and the colour shows no 
significant change during the transit. The measured planet- 
star radius ratios for the two bandpasses are also consistent 
within 3 (t (Table |4|. The FPP computed from the Kepler 
light curve is 1.91%. Given the low FPP from Kepler and 
that we measure no significant difference in the planet-star 
radius ratios, we believe this candidate is likely a planet. 
However, we encourage additional follow-up to confirm these 
findings. We discuss this target further in i)5.1l 

KOI 531.01: Systematics are clearly present in the light 
curve prior to the transit, but we find no correlation between 
the systematic trends seen and either instrumental or astro- 
physical parameters. Despite the systematics, we find that 
the planet-star radius ratios differ in the two bandpasses by 
>3a (Table [4|, suggesting this to be a false positive. Also, 
the FPP we compute from the Kepler light curve is very 
high (90.0%). However, we only measure a 1.5cr difference in 
the colour during transit compared to out-of-transit, which 
makes the identity of this target somewhat ambiguous. We 
discuss these results further in 

KOI 732.01: The light curve suffers from low signal-to- 
noise due to the faintness of the target (Kp = 15.3), but 
we find no visual evidence of a wavelength-dependent tran¬ 
sit depth. However, the measured planet-star radius ratios 
differ by >3a (Table|4]). We note that the scatter in the color 
measured for KOI 732 is significantly smaller than the scat¬ 
ter for our other targets. This could be due to particularly 
favorable observing conditions and/or unusually low stellar 
activity for KOI 732. Given the relatively small number of 
observations used to estimate the photometric precision, we 
are cautious not to over-interpret these results. We recog¬ 
nize that the scatter of photometry might be smaller than 

® While the depth in the Kepler bandpass is shallower than in 
the GTC bandpasses, we attribute this to dilution resulting from 
a nearby star that contaminated the Kepler aperture. 
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Figure 1. Transit light curves, residuals, and color for the 2012 
June 11 UT transit of KOI 439.01. The blue and red points repre¬ 
sent the 666 nm and 858 nm data, respectively. The solid curves 
are the best-fit models. Representative photometric error bars are 
shown on the upper left-hand side of the top and bottom panels. 
The dashed line indicates the depth of the transit measured in the 
Kepler bandpass (Table (2]). The dash-dot line in the top panel 
illustrates the quality of the light curve model. The dash-dot line 
in the bottom panel indicates where the color is zero (for refer¬ 
ence). The errors on the color are significantly smaller than the 
individual photometric errors shown in the top panel as a result 
of binning the data prior to computing the color. 
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Figure 2. Same as Figure [T] but for the 2012 September 5 UT 
transit of KOI 531.01. 


expected given the relatively small number of observations. 
Furthermore, the low FPP of 0.544% measured from the 
Kepler data support our visual findings that there is no sig- 
nihcant depth difference between our bandpasses. Overall, 
we conclude that this is likely a real planet (which is further 
supported by the multiplicity of this system; 



- 3 - 2-10 1 2 3 
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Figure 3. Same as Figure[T] but for the 2012 May 27 UT transit 
of KOI 732.01. 


Table 4. Median Ratios of Planet-Star Radius Ratios 


KOI 

P666/P858 

'^peee/psss 

439.01 

0.975 

0.082 

531.01 

0.928 

0.046 

732.01 

0.943 

0.004 


“ The 3(7 value is presented here. 


5 DISCUSSION 
5.1 KOI 439.01 

iLillo-Box et al.l (l2012l f identihed a companion 5.34 mag 
fainter in i-band and 5.453 arcsec from the target. As our 
aperture was 3. 0-3.3 arcsec , we a voided contamination from 
this companion. iLaw et al. IlM) also ruled out companions 
within 2.5 arcsec of the target. From this, we believe that 
the transit signal is indeed coming from the target. Our FPP 
calculations support this as well (Tableland However, 
we recommend radial velocity as well as additional photo¬ 
metric follow-up to help confirm the nature of this target. 
Such observations could allow a background eclipsing binary 
scenario to be confidently ruled out. 

It should be noted that KOI 439.01 was included in the 
third version of the Kepler eclipsing binary catalog^ Fur¬ 
thermore, a second planet candidate was found in this sys¬ 
tem with a period of 5.4 days and a depth of 117 ppm (KOI 
439.02), but was later fla gged as a false alarm (iLissauer et al.l 
|2014 iRowe et al.l [ 2 ^ . While unrelated, these facts do 
raise a cautionary flag towards the planetary nature of KOI 
439.01. Regardless, our observations show no obvious evi¬ 
dence that this target is an eclipsing binary false positive or 
false alarm, strengthening the case that KOI 439.01 is likely 
a planet. _ 

This target was noted bv lBeauee fc Nesvornvl ll2013l) as 


® However, it is our understanding that the third version of 
the eclipsing binary catalog includes all possible eclipsing binary 
stars, and it is not a comprehensive catalog of confirmed Kepler 
eclipsing binaries; http://keplerebs.villanova.edu/ 
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Table 3. Best-Fit Model Parameters 


KOI 

*0 

(BJD_TDB-2455000) 

T 

(d) 

b 

Rp/Ri, 

(666 nm) 

Rp/R* 

(858 nm) 

Cl 

(666 nm) 

C2 

(666 nm) 

Cl 

(858 nm) 

(S 

439.01 

1089.60915±2.0796e-04 

0.0912±1.30e-03 

0.0000±7.231e-02 

0.0513±1.18e-03 

0.0526±8.45e-05 

0.7099 

0.3153 

0.5536 

( 

531.01 

1176.49398±8.6418e-06 

0.0413±1.29e-04 

0.9675±2.133e-03 

0.0809±1.01e-03 

0.0872±9.33e-04 

0.7212 

-0.0038 

0.6056 

- 

732.01 

1074.62127±3.8018e-05 

0.0630±1.69e-04 

0.0000±4.507e-03 

0.0300±2.25e-05 

0.0318±3.95e-05 

0.6887 

0.1917 

0.5583 

- 


Table 5. False Positive Probabilities 


KOI 

Ac (nm) 

PpL 

Peb 

Pheb 

Pbeb 

FPP 

439.01 

Kepler 

0.981 

2.94e-5 

5.56e-13 

0.0191 

0.0191 

531.01 

Kepler 

0.100 

0.857 

1.87e-4 

0.0428 

0.900 

732.01 

Kepler 

0.995 

2.72e-15 

6.33e-34 

5.44e-3 

5.44e-3 


a particularly interesting planet, if confirmed. This is be¬ 
cause based on the estimated planet radius (3.9 i?®) and 
the short orbital period (1.9 days), this candidate sits within 
what is becoming known as the sub-Jovian desert. We dis¬ 
cuss this region in further detail in >15.51 


5.2 KOI 531.01 


This candidate has been of particular interest to groups 
studying planets around cool stars, as the estimated stel¬ 
lar effective temperature is •^4000 K llMuirhead et al.l 


20121: Dressing & i Charbonneaull2013l : lMorton fc Swift 20 14 


Swjft_etj^ 20jjl) . Wh ile commonly treated as a planet. 


Morton fc Swift ll20l4l notably find KOI 531.01 to have 


a calculated false positive probability of 99%. A more re¬ 
cent calculation using all available quarters of Kepler data 
yield ed a false positive probability of 48% ll Swift et al.l 
l2015h . which is still high but not conclusive evidence of a 
false positive nature . Furthermore, this target was noted in 
iBorucki et al.l (l201lh as having a “strange” light curve. The 
phase-folded Kepler light curve presented in the data vali¬ 
dation summary from the NASA Exoplanet Archive shows a 
clear and fairly v-shaped transit. The “strangeness” appears 
to come from several points during transit that appear to be 
mis-phased and possibly suggest the presence of a secondary 
eclipse, as they are shallower than other points during tran¬ 
sit. However, the latest light curve generated using data from 
Q1-Q17 shows no outliers and no sign of a secondary event. 
Thus, the only published observatio nal evidence of thi s can¬ 
didate being a false positive is from I Swift et al.l (l2013ll . who 
claim to have found a deep secondary eclipse for KOI 531.01 
in the Kepler data (which appears to contradict the latest 
light curve from the Kepler team). While the raw Kepler 
light curve may hint at a potential depth difference in odd 
vs. even transits, it is not definitive. Furthermore, while KOI 
531.01 is also included in the latest Kepler eclipsing binary 
catalog, the light curve in the catalog is not particularly 
informative and does not display evidence of a secondary 
eclipse or depth differences. La stly, some evidenc e of transit¬ 
timing variations was found bv iFord et al.l (1201 ll L but it was 
noted that the short duration likely affected the transit times 
so the variations may not be significant. 

The depth difference measured in our photometry (Ta¬ 


ble |3J co mbined with the high FPP computed by our group 
(and by iMorton fc SwiftI 1201 J and ISwift et al.ll2015f) from 
the Kepler photometry suggests that KOI 531.01 is indeed 
a false positive. Our FPP calculations specifically suggest 
that KOI 531.01 is itself an EB (rather than being a hi¬ 
erarchical system or a blend with a background eclipsing 
system). This is supported by iLaw et ^ I 2 OI 4 I I , who de¬ 
tected no c ompanions down to Am = 5 within 2.5 arcsec. 
Given that iLaw et all ll2014h was sensitive only to targets 
up to ~5 mag fainter at distances of 0.5 arcsec from the 
target, it is still possible that there is a faint stellar com¬ 
panion blended with the host star that either affected our 
photometry and/or is the source of the observed transit sig¬ 
nal. Ultimately, because the measured color difference has 
only a marginal significance (1.5 ct), we declare that the na¬ 
ture of KOI 531.01 remains somewhat ambiguous. We urge 
radial velocity and additional photometric follow-up (espe¬ 
cially at redder wavelengths) to help clarify the nature of 
this system. If real, it would be one of the few known giant 
(~6 I?©) planets with an M dwarf host. 


5.3 KOI 732.01 

After our observations were conducted, the Kepler team 
identified two additional candidate planets around KOI 
732. Given the extremely low false p ositive probabil i ty for 
candidates in m ulti-planet systems (iLissauer et al.l l2014l : 
iRowe et al.ll20l3 f. this planet has already effectively been 
validated. Our findings support this argument, as does the 
low FPP computed from the Kepler light curve (Table [5]). 

We note that the phase-folded Kepler light curve given 
in the data validation summary from the NASA Exoplanet 
Archive for this target displays ‘brightening’ features before 
ingress and after egress. We find no evidence of this in our 
light curve, and we attribute this to be an artifact of the 
reduction process in the Kepler pipeline. 

This is another target that was noted by 
iBeauee fc Nesvornvl ll2013ll as being within the sub-Jovian 
desert, which we discuss further in >15.51 
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5.4 Kepler False Positive Rate 


In this paper we present observations that support the case 
for a planetary nature for KOI 439.01 and KOI 732.01 
and a likely false positive nature for KOI 531.01. Follow- 
iuR the hiah false posit ive p robabilities presen ted here and in 
iMorton fc SwiftI (120141 ') and lSwift et al.f ((201^ , we treat KOI 
531.01 as a false positive in the discussion below. In Papers 
I and II, we identified a total of three false positives (KOI 
565.01, 225.01, and 1187.01) and supported two candidates 
as validated planets (KOI 420.01 and 526.01). Thus, of eight 
candidates observed in total, we identify four as viable plan¬ 
ets and identify four as false positives. This suggests a false 
positive rate as high as 50% for short-period candidates. We 
note that KOI 565.01 and KOI 1187.01 were identified as 
false positives upon further analysis by the Kepler team, so 
if we exclude these from our sample as well as KOI 531.01 
we find that one of five or 20% of short-period candidates 
are false positives. While our sample size is small, the higher 
false positive rate w e estimate is consistent with results from 
iFressin et al.l (l2013l i for short-period candidates. Lower false 


sample of candidates (e.g. 

Morton & Johnsonll201ll: iMortonI 

l20li IFressin et al.l l2013l: 

Desert et al.| 2015|). while others 


have estimated false positi ve rates as high as ^3 5% for giant 
planets with P < 25 days (ISanterne et al.ll2 012 ) or for can¬ 
didat es due to some type of contamination ( Coughlin et al.l 


In FigurelH we present the latest distribution of Kepler 
planet candidates, eclipsing binaries, and false positives^ 
The false positives are those KOIs that are listed with a 
Kepler disposition of FALSE POSITIVE. We note that 
the cumulative catalog was used to generate these distri¬ 
butions, which has not been fully vetted to date. Therefore, 
some false positives could exist among the KOIs listed as 
candidates. Notably, KOI 225.01 is currently listed in the 
archive as a planet candidate, but it is clearly a false pos¬ 
itive based on the significant wavelength-dependent transit 
depth (Paper II). KOI 531.01 is also listed as a planet can- 
didate, though it has a false posit i ve pr obability 50% 
jMorton &: Swifj l2014l : ISwift et al.l l2015h . These distribu¬ 
tions clearly emphasize the continued prevalence of eclipsing 
binaries (and correspondingly, false positives) at the short¬ 
est orbital periods, as discussed in further detail in Paper 


11 . 


Given the small sample size presented here, we cannot 
make strong statements regarding trends in the false positive 
rate with other planetary or stellar parameters. Still, after 
investigating different parameters that included the Kepler 
magnitude and Galactic latitude, we find no evidence of cor¬ 
relations in the properties of our identified false positives or 
their hosts. However, we find that at short periods {P ^ 6 
days), the distribution of KOIs dispositioned as false posi¬ 
tives is clearly skewed towards lower Galactic latitudes com¬ 
pared to the candidate population (Figure O. This is to be 
expected, given the increased density of stars closer to the 
Galactic plane. Somewhat surprisingly, the distributions of 


^ Accessed from the NASA Exoplanet Archive 
(http://exoplanetarchive.ipac.caltech.edu/) on 2015 

April 4 and the Kepler eclipsing binary catalog 
(http://keplerebs.villanova.edu/) on 2014 June 25. 



Figure 4. Distribution of Kepler planet candidates (green), false 
positives (red), and eclipsing binaries (blue) as a function of or¬ 
bital period. The arrows indicate the periods of the eight KOIs in 
our entire sample. The green arrows indicate KOIs we tentatively 
identify as planets, and the red arrows indicate those KOIs iden¬ 
tified as false positives. The two red arrows filled with green are 
KOIs that are still listed in the Kepler catalog as active planet 
candidates (KOI 225.01 and KOI 531.01, though the false positive 
nature of the latter is not definitive). 



Figure 5. Similar to Figure|4l but for the distribution of Kepler 
planet candidates (green) and false positives (red) as a function 
of Galactic latitude. Only planet candidates and false positives 
with periods less than 6.5 days are included. The arrows are the 
same as in Figure |3] 


Kepler-identified false positives and candidates in terms of 
their Kepler magnitude are similar, except at the faintest 
magnitudes (Kp > 16) where false positives are more preva¬ 
lent (also as expected). 

Lastly, we do not expect our findings to be otherwise 
biased by our selection criteria, though we note that our re¬ 
striction on transit depth and/or duration could technically 
bias the spectral types we target. We do not expect such 
a bias in our current target sample, given that the stellar 
effective temperatures of our targets range from about 4000 
to 6000 K. 
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5.5 Planets in the sub-Jovian Desert 


With the increasing numbers of confirmed exoplanets and 
the thousands of planet candidates from Kepler, exciting 
trends have been emerging in the planet population. In par¬ 
ticular, there is a region in the period-radius plane (and cor¬ 
respondingly in the period-mass plane) where there appears 
to be no known planets. This region has been referred to as 
the sub-Jovia n desert and is the sub j ect of several recent 
studi e s (e.g., Bem'tez-Llambav et al. 201ll;_ Szabo_&J[^ 


20n; iBeauge fc Nesvornvl 1201 iKurokawa fc Nakamotc 
2014 1. The region encompasses orbital periods approxi¬ 
mately less than 2.5 days and radii between about 3 and 11 
(Figure ini). For comparison purposes, we show in Fig¬ 
ure [7] the desert in the period-mass plane (where it spans 
approximately 0.03—1 Mj). 


We note that in Figs [6] and [71 the desert region is 
overemphasized because we show planets discovered by both 
radial velocity surveys and ground- and space-based transit 
surveys. This combined sample therefore suffers from many 
biases. Still, the lack of planets in this region cannot be 
easily explained by observational bias. Short-period planets 
are generally easier to detect, especially if they are larger 
or more massive than Earth. Assuming the desert is not 
due to observational bias, this suggests that there is some 
formation mechanism at play. For instance, if orbital mi¬ 
gration is required for planets to end up on close-in orbits, 
then there might be a limiting mass where planets can no 
longer migrate inward. The pile-up of hot Jupiters around 
~3 days is suggestive of this. However, since smaller/less 
massive planets are (apparently) allowed to orbit very close 
to their host stars (e.g. with periods less than 1 day), it 
has been theorized that these small planets are the leftover 
cores of planets th at had their a tmospheres blown away by 
the host star (e.g. lYoudinll2011^ . This mechanism is likely 
most effective for planets with low surface g r avity ^ and not 
necessarily for all hot Jupiters. I Owen fc Wu (|2013h investi¬ 
gated evaporation in detail and conclude that evaporation 
explains a lack of low density planets on close-in orbits as 
well as a deficit of planets around 2 R^. Currently, this 
seems to be the most viable explanation to explain the ob¬ 
served distribution of pl anets and candidates. While beyond 
the scope of this paper, iBeauge fc Nesvornvl (|2013l ') offer a 
more complete discussion of alternative causes of the desert, 
including that hot Jupiters may have been tidally captured 
while super-Earth-size planets may have been affected by 
disk-planet interactions. 


Given the competing formation mechanisms, any can¬ 
didate planets within this desert are therefore of extreme 
interest. In total, there are currently ~30 KOIs within the 
desert, including KOI 439.01 and KOI 732.01 which we iden¬ 
tify as likely planets here. However, several KOIs have prop¬ 
erties that place them near the boundaries of the region, and 
uncertainties in the planetary radius can easily shift some 
of them out of the desert. Follow-up observations of these 
KOIs are needed to confirm their planetary nature. 

At least three of the KOIs in the desert are in known 
multiple candidate systems, suggesting that these are likely 
real planets. Along with our observations providing addi¬ 
tional evidence in support of KOI 439.01 and 732.01 as 
viable planets, and three other confirmed planets (Kepler- 



Orbitol Period [days] 


Figure 6. Planet radius versus orbital period for all confirmed 
transiting planets (from both ground-based surveys and Kepler-, 
gray squares) and for active Kepler planet candidates (green cir¬ 
cles). The black squares mark the locations of the four KOIs 
we support as planets, while the red circles mark the four KOIs 
that we identified as false positives. We consider KOI 439.01 and 
732.01 to be planets and KOI 531.01 to be a false positive in 
this context. The parameters for confirmed planets were taken 
from The Exoplanet Orbit Database at exoplanets.org on 2014 
July 31. The region marked by dashed black lines indicates the 
regime of the so-called sub-Jovian desert (for Rp ~ 3—11 and 
P < 2.5 days). While we initially selected KOIs with radii les s 
than 6 i?®, updated stellar parameters from iHuber et ahl (l20l4 f 
and light curve models from the Kepler team yielded new planet 
radii of 33.63 and 31.69 R^ for the false positives KOI 1187.01 
and 225.01 (the two shortest period KOIs we observed, the lat¬ 
ter of which is notably still included in the archive as a planet 
candidate). 


41b, Kepler-119b, WASP-43b)|3 there could be at least eight 
known planets located in the desert so far (albeit towards 
the edges of the region). It will be interesting to see if the 
KOIs that are in the more central region of the desert sur¬ 
vive the vetting process, considering the higher false positive 
rate that exists for KOIs with short periods compared to the 
global rate. It may be that planets can exist in this regime, 
but under atypical circumstances of planetary formation and 
migration. If none do, it at least appears that some planets 
are filling out the edges of the region and as a result, the 
boundaries of the desert are being constrained. The con¬ 
tinued study of planets in this region is therefore critical for 
informing planetary formation theories and at the very least 
will help to constrain the boundaries of the desert. 


6 CONCLUSION 

In this paper, we use multi-wavelength transit photometry 
to argue for a planetary nature for the two Kepler planet 
candidates KOI 439.01 and KOI 732.01 and a false positive 
identification for KOI 531.01. Additional follow-up observa¬ 
tions, such as radial velocity and/or photometry at redder 
wavelengths, are needed to clarify the nature of KOI 439.01 


® Determined from parameters given on exoplanets.org 
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Figure 7. Similar to Figure|6l but showing the minimum planet 
mass versus orbital period for all confirmed planets (from both 
radial velocity and transit surveys). The region marked by dashed 
black lines covers the mass range of 0.03—1 Mj and P < 2.5 days. 


in particular. In combination with r esults from previo usly 
published work llColon &: FordlEoill : IColon et al.l]2012l) . we 
have vetted a total of eight small, short-period Kepler can¬ 
didates through GTC observations. While our sample size is 
small, our results support a slightly higher false positive rate 
for short-period planet candidates than for the overall sam¬ 
ple. This is in agreement with results from other studies and 
has interesting implications for Kepler planet candidates lo¬ 
cated within the sub-Jovian desert, a region in the period- 
radius plane where some planetary formation mechanisms 
appear to restrict planets from living there. Future obser¬ 
vations of these candidates will improve our understanding 
of different formation scenarios, particularly if none of these 
candidates survive the vetting process. 
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